Magnetotransport of Dirac Fermions on the surface of a topological insulator 
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We study the properties of Dirac fermions on the surface of a topological insulator in the presence 
of crossed electric and magnetic fields. We provide an exact solution to this problem and demonstrate 
that, in contrast to their counterparts in graphene, these Dirac fermions allow relative tuning of 
the orbital and Zeeman effects of an applied magnetic field by a crossed electric field along the 
surface. We also elaborate and extend our earlier results on normal metal-magnetic film-normal 
metal (NMN) and normal metal-barrier-magnetic film (NBM) junctions of topological insulators 
[Phys. Rev. Lett. 104, 046403 (2010)]. For NMN junctions, we show that for Dirac fermions with 
Fermi velocity vf, the transport can be controlled using the exchange field J of a ferromagnetic 
film over a region of width d. The conductance of such a junction changes from oscillatory to a 
monotonically decreasing function of d beyond a critical J which leads to the possible realization 
of magnetic switches using these junctions. For NBM junctions with a potential barrier of width d 
and potential Vo , we find that beyond a critical J , the criteria of conductance maxima changes from 
X = eVod/hvF = nil to \ — ( n + 1/2) 71 " for integer n. Finally, we compute the subgap tunneling 
conductance of a normal metal-magnetic film-superconductor (NMS) junctions on the surface of a 
topological insulator and show that the position of the peaks of the zero-bias tunneling conductance 
can be tuned using the magnetization of the ferromagnetic film. We point out that these phenomena 
have no analogs in either conventional two-dimensional materials or Dirac electrons in graphene and 
suggest experiments to test our theory. 

PACS numbers: 71.10.Pm, 73.20.-r 
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I. INTRODUCTION 

Topological insulators with time reversal symmetry in 
two and three dimensions (2D and 3D) have been stud- 
ied extensively in recent years, both theoretically and 
experimentally^ - — . The 3D topological insulators can be 
characterized by four integers uq and ^1,2,3^- The first in- 
teger specifies the class of topological insulators as strong 
(^0 = 1) or weak (i/q = 0), while the last three integers 
characterize the time-reversal invariant momenta of the 
system given by Mq — (v%b\, vibz)/2, where 61,2,3 
are reciprocal lattice vectors. The topological features 
of strong topological insulators (STI) are robust against 
the presence of time- reversal invariant perturbations such 
as disorder and lattice imperfections. It has been theo- 
retically predicted^ and experimentally verified 2 that 
the surface of a STI has an odd number of Dirac cones 
whose positions are determined by the projection of Mo 
on to the Brillouin zone of the surface. The position and 
number of these cones depend on both the nature of the 
surface concerned and the integers ^1,2,3- For compounds 
such as HgTe and Bi2Sc3, specific surfaces with a single 
Dirac cone near the T point of the 2D Brillouin zone have 
been found2i££. Such a Dirac cone is described by the 
Hamiltonian 



H 



i)\k) (hv F a ■ k - ip[k), (1) 



where <r(i) denotes the Pauli (identity) matrices in spin 
space, ip = (^,ipi) T is the annihilation operator for the 
Dirac spinor (T denotes the transpose of a row vector), 



vp is the Fermi velocity, and /i is the chemical potential 9 . 
Recently, several novel features of these surface Dirac 
electrons such as the existence of Majorana fermions in 
the presence of a magnet-superconductor interface on 
the surface&i^, generation of time-reversal symmetric 
Px + ipy-wave superconducting state via proximity to a s- 
wave superconductor—, anomalous magnetoresistance of 
ferromagnet-ferromagnet junctionsAi', and novel spin tex- 
tures with chiral properties^ have been studied in detail. 
Further it has been shown in Ref. [ID that it is possible 
to realize a magnetic switch by magnetically tuning the 
transport of Dirac fermions with a proximate ferromag- 
netic film. However, the response of these fermions in 
the presence of crossed electric and magnetic fields has 
not been studied so far. Another aspect of such fermions, 
namely, their transport through a normal metal-magnetic 
film-superconductor (NMS) junction has also not been 
explored. 

In this work, we study several magnetotransport prop- 
erties of these surface Dirac fermions in experimentally 
realizable situations. We first study the properties of 
the fermions in the presence of crossed magnetic [B = 
(0, BcosO, BsinO)] and electric fields [£ — (5,0,0)] as 
shown in Fig. [TJ We present an exact solution of this 
problem and show that for j3 = £ / (vfB smO) < 1, 
the relative contributions of the Zeeman and the or- 
bital terms to the Landau level energies, and hence their 
magnetic field dependence, can be tuned by varying ei- 
ther the strength of the applied electric field or the tilt 
of the applied magnetic field. We also show that for 
£ > vfB sin 8, the conductance of these Dirac fermions 
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FIG. 1: Schematic representation of the crossed electric and 
magnetic field geometry. The electric field is applied along x 
while the magnetic field is in the y — z plane. See text for 
details. 



has an unconventional dependence on the tilt angle 9 of 
the applied magnetic field and that this dependence can 
be used to realize electric-field controlled switching. 

The second study involves an extension of the re- 
sults obtained in Ref. [l2| regarding transport in normal 
metal-magnetic film-normal metal (NMN) and normal 
metal-barrier-magnetic film (NBM) junctions of topolog- 
ical insulators. The relevant experimental geometries are 
shown in Fig. [5] We study the transport of these Dirac 
electrons across a region with a width d where there is 
a proximity-induced exchange field J arising from the 
magnetization fh — nioy of a proximate ferromagnetic 
film as shown in the left panel of Fig. [2] We demonstrate 
that the tunneling conductance G of these Dirac fermions 
through such a junction can either be an oscillatory or 
a monotonically decaying function of the junction width 
d. One can interpolate between these two qualitatively 
different behaviors of G by changing mo (and thus J~) 
by an applied in-plane magnetic field leading to the pos- 
sible use of this junction as a magnetic switch. We also 
study the transport properties of Dirac fermions across 
a barrier characterized by a width d and a potential Vq 
in region II with a magnetic film proximate to region 
III as shown in the right panel of Fig. [5J We note that 
it is well known from the context of Dirac fermions in 
graphenei^ that such a junction, in the absence of the 
induced magnetization, exhibits transmission resonances 
with maxima of transmission at x = eVod/hvF = nit, 
where n is an integer. Here we show that beyond a criti- 
cal strength of mo, the maxima of the transmission shifts 
to x = (n + l/2)7r. Upon further increasing mo, one can 
reach a regime where the conductance across the junc- 
tions vanishes. We also point out that such NMN and 
NBM junctions can be used to determine the exact form 
of the Dirac Hamiltonian on the surface of the topological 
insulator. 

Finally, we study the transport of Dirac fermions 
across a NMS junction as shown in Fig. [3J The interme- 
diate region (region II) in this junction has a thickness 
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FIG. 2: Proposed experimental setups. Left panel: Schematic 
representation of a NMN junction. The ferromagnetic film 
extends over region II of width d providing an exchange field 
in this region. Right panel: Schematic representation of a 
NBM junction. The film extends over region III while the 
region II has a barrier characterized by a voltage Vq. V and 
I denote the bias voltage and current across the junctions 
respectively. See text for details. 



d with a proximate ferromagnetic film providing a mag- 
netization M, while superconductivity is introduced in 
region III via the proximity effect. We provide a detailed 
analysis of the subgap tunneling conductance G of these 
NMS junctions as a function of the applied voltage V 
and magnetization M. In particular, we point out that 
the positions of the maxima of the zero-bias tunneling 
conductance in such NMS junctions as a function of the 
width d of the magnetic film can be varied by tuning the 
induced magnetization M. We stress that the proper- 
ties of the Dirac fermions elucidated in all these studies 
are a consequence of their spinor structure in physical 
spin space, and thus have no analogs for either conven- 
tional Schrodinger electrons in 2D or Dirac electrons in 
graphenei^— . 

The organization of the rest of the paper is as follows. 
In SecHH we study the properties of the Dirac fermions in 
the presence of crossed electric and magnetic fields. This 
is followed by the study of NMN and NBM junctions of 
these Dirac materials in Sec. lIIII In Sec. lIVl we study the 
transport properties and subgap tunneling conductance 
of NMS junctions of topological insulators. Finally we 
discuss possible experimental verification of theory and 
conclude in Sec. W\ 



II. CROSSED ELECTRIC AND MAGNETIC 
FIELDS 

We begin with the properties of Dirac electrons in a 
crossed electric and magnetic field as shown in Fig. [T] 
The Hamiltonian for the Dirac Fermions for this case 
can be written as 

H = J d 2 r^(f)[vF<J ■ n — [il — gnB& ■ B — e£ x]ip(r), 

(2) 

where n = — ihV — eA is the canonical momentum, c 
is set to unity, g is the gyromagnetic ratio, fig is the 
Bohr magneton, and we choose the vector potential to be 
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FIG. 3: Schematic representation of the NMS junction. The 
ferromagnetic film is on region II while superconductivity is 
induced in region III via a proximate superconducting film. 
See text for details. 



A = (0, Bxs'md, — Bxcosd). Note that here the Zeeman 
term does not determine the spin quantization axis of 
the Dirac electrons due to the presence of the II term. 
Thus the in-plane component of the magnetic field, which 
enters the Hamiltonian only through the Zeeman term, 
only provides a constant shift to k y which can be gauged 
away. This property of the Dirac fermions is distinct 
from their counterpart in graphene. For £ = 0, Eq. ^) 
admits a straightforward solution and yields the Landau 
level spectrum 



E n — ±hv F l B 1 V\ n \ + a B sin 6 if n ^ 
= — \gfiBBsm9\ if n = 0, 



(3) 



where a — g 2 [i B /hv F e, Is = y/h/[eB sin 8} is the mag- 
netic length and we also define l B = Ib{9 = tt/2) for 
later use. The n — state is non-degenerate as is 
also known from analogous studies of Landau levels in 
graphene^. For the Dirac electrons on the surface of 
HgTe, v F ~ 5.5 x 10 5 m/s, so that a ~ 1(T 4 /T leading 
to a negligible contribution of the Zeeman term in the 
spectrum. 

The situation changes when an electric field is applied 
along x. In this case, for £ < v F B sin (9, one can define a 
boost parameter j3 = £ /[v F B sind] < 1 and carry out a 
Lorentz transformation^ 

x = x, y' = 7(2/ - pv F t), t' = j(t- py/v F ), 

£' = j(£-f3B sinO), 

B' sin & = 7(5 sin 6>- (3£), 

tp'(f') = exp[-er y arctanh(/3)/2] i/j(r), (4) 

where 7 = (1 — /3 2 )" 1 / 2 . In the boosted frame the 
Schrodinger equation reads 



■1 I a 1 fa .eB'sinO' s 
imip I <j x o x + <Jy{o y i — 1 X, 



gUB&zB sinO ip'. 



(5) 




FIG. 4: Plot of Landau level energy E„ for n = I as a function 
of the tilt angle 6 and transverse momentum k y for a fixed 
electric field £ — 0.5vfB. Note that 9 = tt/6 corresponds to 
13 = 1. 



Note that such a boost transformation affects the orbital 
part of the magnetic field only; the Zeeman field remains 
unchanged. The energy eigenvalues of Eq. ([5]) can be 
easily obtained and are given by 



EL 



±Hv F l B 1 y f[r 
-\g^ B Bsin9 



+ aB sin 6 if n ^ 
if n = 0, 



(6) 



where l B = ^/h/ [e£?'sin#']. Then a reverse boost to the 
"laboratory" frame yields 



E n (k y ) 



±hv F i B ln / 



aBj sin( 



3/2 a/W 

Phvpky if n =^ 
^~ x \g[iBB sin#| — (ihv F k y 



ifn = 0.(7) 



Eq. ([7]) is one of the central results of this section. 
It demonstrates that a collapse of the Landau levels for 
the Dirac fermions can be induced by varying either the 
electric field for a fixed tilt of the applied magnetic field 
or by varying the magnetic field tilt for a fixed electric 
field. A plot of the energy level as a function of this tilt 
and the transverse momentum k y is shown in Fig. |4] 

We also note that the magnetic field dependence of 
the Landau level energy gaps are different from their 
counter part in graphene. To illustrate this we define 
A„ = E n+1 (k y ) - E n {k y ) for n ^ 0. For /5 « 1, we find: 



A, 



hv F l° B 



(8) 



For j3 ~ 1 when jaB \n\, we find that 

AJinvplg- 1 ) ~ 7 - 3 / 2 /(2 v / ^t). This behavior 
is distinct from its counter part in graphene where 
= 7 - 3 / 2 VsTri0(V^" 



A n /(hv F l 



1 - y/n) for all B. 
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FIG. 5: Plot of Landau level energy gap A„ for n — 1 (black 
solid line) 2 (red long-dashed line) and 3 (blue dashed line) 
as a function of the tilt angle 9 for a fixed electric field £ = 
0.2v F B 



G/Gn 




e 



0.5 1 1.5 

FIG. 6: Plot of the conductance G/Go as a function of the 
tilt angle 6 for electric field B = 1 (black solid line), 1.2 (green 
dash-dotted line), 1.5 (red long dashed line), 3 (blue dashed 
line) for a fixed magnetic field B. 



However, since a ~ 10~ 4 /Tesla, this behavior can only 
be seen in a very tiny window near critical tilt and would 
be hard to figure out in experiments. On the contrary, 
the variation of dispersion of A„ with the tilt angle 8, 
seen in Fig. [S] for several n, can be tested experimentally 
with microwave absorption experiments routinely done 
for conventional quantum Hall systems^. This will be 
discussed further in Sec. [V] 

Now we turn to the solution of this problem in the 
regime £ > vfB where we get scattering states. In this 
regime, we define a parameter ft' = vfB sm0/£ and per- 
form a similar boost transformation as outlined earlier. 
This allows us to shift to a reference frame where there 
is no magnetic field and the Schrodinger equation, in the 
momentum representation, readsi^ 



e£ Kdk x ip = hvp(a x k x + a y k y ) — E' 
—gfiBPzB sin# ip', 



(9) 



where £' = j'(£ - ft'v F B sin0) and E' = j'(E-hv F k y ft') 
are the electric field and energy as seen in the boosted 
frame, and 7' = 1/ yjl — ft' 2 . The scattering states 
can now be easily obtained from this equation by not- 
ing the similarity of this equation with the standard 
Landau-Zenner problem with modified Planck's constant 
h — > e£'h. In particular, the transmission probability 
of these Dirac electrons in the direction of the applied 
electric field in the boosted frame can be written as 



T(k y ;E) 



-irdl~/'[( gf j, B Bsin8) 2 + (hv F k y ) 2 ]/(hv F ) 2 



-(10) 



where d a = y/ Hvf / (e£) is the legth scale set by the elec- 
tric field. Now one can rotate back to the "laboratory 
frame" using k y = j'iky — ft'e/HvF), and integrate over 
k y modes to obtain the tunneling conductance^ 



G 



G (l-ft 2 ) 



' 2)3/4 e -ny'dl{giJ.B B sin 8) 2 /(hv F ) 2 



(11) 



where Go = e 2 L y /(hdo), and L y is the sample width. 
We find that in contrast to graphenei^, the Zeeman term 
arising from the magnetic field along z produces an addi- 
tional exponential suppression of the conductance. This 
can be understood by noting that in topological insu- 
lators, a Zeeman magnetic field along z results in the 
generation of a mass term for the Dirac electrons and 
hence leads to a suppression of the conductance. A plot 
of G/Gq as a function of the tilt of a magnetic field 6 is 
shown in Fig. [6] for several representative values of the 
electric field £ and for a fixed magnetic field B. The plot 
shows that for small electric fields, the conductance is 
quickly suppressed as we increase from to 7r/2; how- 
ever for larger fields, the suppression is minimal. Thus 
one can tune the conductance of these insulators either 
by tuning the electric field at a fixed 8 or by tuning 9 at 
a fixed electric field. 



III. TRANSPORT IN NMN AND NBM 
JUNCTIONS 

In this section, we analyze the properties of NMN and 
NBM junctions of topological insulators as shown in the 
left and right panels of Fig. [2l Sec. IIII Al discusses the 
NMN junctions while Sec. IIII Bl elucidates the properties 
of the NBM junctions. 



A. NMN junctions 

The proposed experimental set up for the NMN junc- 
tion is shown in the left panel of Fig. [2j The Dirac 
fermions in region I and III are described by the Hamil- 
tonian in Eq. (JlJ. Consequently, the wave functions of 
these fermions moving along ±2; in these regions for a 
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fixed transverse momentum k y and energy e can be writ- 
ten as 



1.0 



ipf = (l,±e ±ia ) exp[i{±k x x + k y y)}, 
where j takes values I and III, and 

a = a,rcsm(hv F k y /\e + ji\), 



(12) 



k x {e) = ^[{e + Li)lhv F Y~kl. 



(13) 



In region II, the presence of the ferromagnetic strip with 
a magnetization mo = moy leads to the additional term 



induced 



dxdy J9(x)9(d- x)i)\r)a y ip(r), (14) 



where J ~ m is the exchange field due to the pres- 
ence of the strips, and 8(x) denotes the Heaviside step 
function. Note that -ffinduced may be thought as a vector 
potential term arising due to a fictitious magnetic field 
Bf = (J /ev F )[5(x) — 5(d — x)]z. This analogy shows 
that our choice of the in-pane magnetization along y is 
completely general; all gauge invariant quantities such 
as the transmission probability are independent of the 
x-component of m in the present geometry. We empha- 
size that this effect is distinct from that due to a finite 
z component of mo which provides a mass to the Dirac 
electrons. For a given mo, the precise magnitude of J 
depends on the exchange coupling of the film and can be 
tuned, for soft ferromagnetic films, by an applied field 11 . 
The wave function for the Dirac fermions in region II 
moving along ±x in the presence of such an exchange 
field is given by 



= (l.ie*'") exp^ifc^ + ^y)], 



(15) 



where 



/3 = arcsin(ftw F (fc y + M)/|e + fx\), M = J/(hv F ), 



= J[{e + ^)/hv F Y-{k y + Mf 



(16) 



Note that beyond a critical M c = ±2|e + /.i\/(Hv F ), and 
hence a critical J c — ±2|e + [t\, k x becomes imaginary 
for all k y leading to spatially decaying modes in region 
II. 

Let us now consider an electron incident on region II 
from the left with a transverse momentum k y and en- 
ergy e. Taking into account reflection and transmission 
processes at x — and x — d, the wave function of the 
electron can be written as 



V-7 



ijjf + rip j , ip. 



ii 



■■ptpu + Qtpin tp. 



in 



III- 
(17) 



Here r and t are the reflection and transmission ampli- 
tudes, and p (q) denotes the amplitude of right (left) 
moving electrons in region II. Matching boundary condi- 
tions on ipj and tpn at x — and tpu and ipm at x = d 
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FIG. 7: Plot of tunneling conductance G/Go for a fixed V 
and fi as a function of the effective width z = d\eV + fi\/hi)F 
for hv F M/\eV + (m\ = 0.3 (green dotted line), 0.7 (blue solid 
line), 1.3 (black dash-dotted line) and 2.1 (red dashed line). 
The value of the critical M is given by hv F M /\eV + fi\ — 2. 



leads to 



te 



1 + r = p + q 

ik x d 



e la - rer la = pe ip 



qe 



-i/3 



_ ik d i ^ — ikd 

pe * + qe * , 
pe i(Kd+l3) _ qe -i{k' x d+p) _ 



(18) 



Solving for t from Eq. (|18)) . one finally obtains the 
conductance 

,tt/2 

G = dl/dV = (G /2) / da T cos a. (19) 

J -7T/2 

Here G = P (eV)we 2 /(irh 2 v F ), p(eV) = \(fi + 
eV)\/[2n(hv F ) 2 } is the density of states (DOS) of the 
Dirac fermions and is a constant for /i 3> eV, w is the 
sample width, and the transmission T = \t\ 2 is given by 

T = cos 2 (a) cos 2 {0)/ '[cos 2 (k' x d) cos 2 (a) cos 2 {fi) 

+ sin 2 (k' x d)(l - sin a sin /3) 2 ]. (20) 

Eq. ([2U]) and the expression for G represent one of the 
main results of this section. We note that for a given 
a, T has an oscillatory (monotonically decaying) depen- 
dence on d provided k x is real (imaginary). Since k x de- 
pends, for a given a, on M, we find that one can switch 
from an oscillatory to a monotonically decaying d depen- 
dence of transmission in a given channel (labeled by k y 
or equivalently, a) by turning on a magnetic field which 
controls mo and hence M. Also, since — 1 < sin a < 1, 
we find that beyond a critical M = M c , the transmis- 
sion in all the channels exhibits a monotonically decay- 
ing dependence on d. Consequently, for a thick enough 
junction one can tune G at fixed V and fi from a fi- 
nite value to nearly zero by tuning M {i.e., mo) through 
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FIG. 8: Plot G/Go versus eV/ fi for several representatives 
values hvFM/n ranging from 3 (left-most black solid curve) 
to 5 (right-most magenta dash-double dotted line) in steps of 
0.5 The effective junction width zo = 5 for all plots. The inset 
shows a plot of eVc/jj, versus hvFM//j,. See text for details. 



M c . Thus such a junction may be used as a magnetic 
switch. These qualitatively different behaviors of the 
junction conductance G for M below and above M c is 
demonstrated in Fig. [7] by plotting G as a function of 
the effective barrier width z = d\eV + h\/Hvf for several 
representative values of fvupM/leV + fi\. Since T and 
hence G depends on M through the dimensionless pa- 
rameter hvpM /\eV + this effect can also be observed 
by varying the applied voltage V for a fixed fi, d, and M. 
In that case, for a reasonably large dimensionless barrier 
thickness zq — d^i/hvp, G/Gq becomes finite only be- 
yond a critical voltage \eV c + fi\ = KvfM/2 as shown in 
Fig- El for several representative values of zq. The critical 
voltage V c can be determined numerically by finding the 
lowest voltage for which G/Go exhibits a monotonic de- 
cay as a function of zq. The plot of eV c //j> as a function 
of HvfM/h, shown in the inset of Fig. [51 demonstrates 
the expected linear relationship between V c and M. We 
note that such a dependence of G on M or V requires the 
Dirac electrons to be spinors in physical spin space, and 
is therefore impossible to achieve in either graphenei^ or 
in a conventional 2D electron gas for which a proximate 
ferromagnetic film would only provide a Zeeman term for 
the electrons, leaving G unaffected. 



B. NBM junctions 

Next, we analyze the NBM junction shown in the right 
panel of Fig. [2] where the region III below a ferromagnetic 
film is separated from region I by a potential barrier in 
region II. Such a barrier can be applied by changing the 
chemical region in region II either by a gate voltage Vq 
or via doping^. We will analyze the problem in the thin 




FIG. 9: Plot of tunneling conductance Gi/Go versus \ an d 
M for a fixed applied voltage V and chemical potential /i. Gi 
vanishes for \M\ > M c = 2\eV + fi\/hv F . 

barrier limit in which Vq — > oo and d — > 0, keeping the di- 
mensionless barrier strength \ = eVod/(hvF) finite. The 
wave function of the Dirac fermions moving along ±x 
with a fixed momentum k y and energy e in this region is 
given by 

rp'jf = (l,± e ±^) exp[i(±Kx + k y y)/V2}, (21) 
where 

7 = &Tcsm(hvFk y /\e + eVo + 
*4'(e) = ^[(e + eVo + ^/hvF} 2 -kl (22) 

The wave functions in region I and III are given by ijij = 
■0/ and ipm = ipiiy where and ipn are given in Eq. 
(|T7| . Note that one can have a propagating solution in 
region III only if |M| < \M C \. 

The transmission problem for such a junction can be 
solved by a procedure similar to the one outlined above 
for the magnetic strip problem. For an electron ap- 
proaching the barrier region from the left, we write down 
the following forms of the wave function in the three re- 
gions I, II and III: ip'j = ijjj +r 1 ipj , ip n = piipjf +q x ijj I J, 
and ipm = ti^n- As outlined earlier, one can then 
match boundary conditions at x = and x = d, and 
obtain the transmission coefficient T\ = \ti\ 2 k x /k x as 

T\ = 2cos/?cosa/[l + cos(/3 — a) 

- cos 2 (x){cos(/3 - a) - cos(/3 + a)}}. (23) 

Note that in the absence of the ferromagnetic film over 
region III, j3 = a, and T x if = cos 2 (a)/[l - 
cos (x)sin 2 (a)]. The expression for T-P, reproduced here 
for the special case of M = 0, is well known from anal- 
ogous studies in the context of graphene, and it exhibits 
both the Klein paradox (Tj 3 = 1 for a = 0) and transmis- 
sion resonances (T° = 1 for x = nw)^. When M ^ 0, we 
find that the transmission for normal incidence [k y = 0) 
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does become independent of the barrier strength, but its 
magnitude deviates from unity: 



normal 



2y/l- [hv F M/\eV + n\Y 
1 + ^/l- (hv F M/\eV + n\) 2 



(24) 



The value of T 1 normal decreases monotonically from 1 for 
M = to for \M\ = \eV + n\/(hv F ), and can thus be 
tuned by changing M (or V) for a fixed V (or M) and \i. 

The conductance of this junction is given by G\ = 
(Go/2) daTiCosa, where are determined from 
the solution of cos /3 = for a given M. A plot of Gi as a 
function of HvpM/\eV + fi\ and x (f° r a fixed eV and /i) 
is shown in Fig. [9j We find that the amplitude of G\ de- 
creases monotonically as a function of \M\ reaching at 
M = M c beyond which there are no propagating modes 
in region III. Also, as we increase M, the conductance 
maxima shifts from \ — n7: to x — ( n + 1/2)tt beyond a 
fixed value of M*(V) ~ icoleT^ + /u|/(/m;.f) as shown in 
the top left panel of Fig. [TUJ Numerically, we find Co = 
0.7075. At M = M*, G x (x = mr) = Gi(x = (n+l/2)7r), 
leading to a period halving of Gi(x) from tt to 7r/2 . 
This is shown in the top right panel of Fig. [TUJ where 
Gi(M = M*) is plotted as a function of x- We note that 
near M*, the amplitude of oscillation of Gi as a function 
of x becomes very small so that G\ is almost independent 
of x- in the bottom left panel of Fig.[TUl we plot x = Xmax 
(the value of x at which the first conductance maxima 
occurs) as a function of hv F M J\eV + fi\ which clearly 
demonstrates the shift. This is further highlighted by 
plotting AGi = Gi(x = 0) — Gi(x = 7r/2) as a function 
of hv F M/\eV + [i\ in the bottom right panel of Fig. [TU] 
AGi crosses zero at M = M* < M c indicating the posi- 
tion of the period halving. Thus the position of the con- 
ductance maxima depends crucially on hv F M/ (eV + fx) 
and can be tuned by changing either M or V. 



C. Alternative forms of the Hamiltonian 

In this subsection, we discuss a possible way of distin- 
guishing between possible forms of the Dirac Hamiltonian 
in the surface of a topological insulator. In the literature 
(see, for instance, Ref. |20I ). two such different forms have 
been studied for the first part of the Hamiltonian in Eq. 
(UJ), namely, 



hi = 
and /12 = 



d 2 k 

d 2 k 
(27T) 



ip^(k) hv F (a x k x + a y ky) ip(k) 
2 ip'{k) hv F (a x k y - (J y k x ) ip(k). 



(25) 



We have implicitly assumed the form hi in the entire 
analysis in this paper. We note that hi and /12 are both 
time-reversal invariant since a — > — a and k — > — k under 
that transformation, and they are also invariant under ro- 
tations in the x — y plane. But under the two-dimensional 




Tiv F M/|eV+n| 



Tiv M/|eV+n| 2 



FIG. 10: Top left panel: Plot of G1/G0 versus x for 
hv F M/\eV + n\ = 0.1 (black solid line), 0.7075 (red dashed 
line), and 1.4 (blue dash-dotted line) for fixed V and /1. Top 
right panel: Plot of G1/G0 versus x at M = M* showing 
the period halving. Bottom left panel: Plot of Xmax versus 
hv F M /\eV + fj,\ showing conductance maxima positions. Bot- 
tom right panel: Plot of AG1/G0 versus hv F M/\eV+fi\ which 
crosses at M = M* . The dotted line is a guide to the eye. 



-x and y — y, they trans- 

~k x ^ ky y ky , 

®x — > &x and 
—hi while /i2 — > Since 



parity transformation x —, 
form differently; since k x 
a y — > —<7 y , we see that hi 
the Hamiltonian of the surface Dirac electrons arises from 
a spin-orbit coupling in the bulk which is then projected 
on to the two-dimensional surface, and we have not dis- 
cussed the bulk Hamiltonian here, we have no a priori 
reason to choose between hi and h^. In principal, we 
could even consider a linear combination of the two such 
as cos# hi +sin# /i2- Clearly, when an in- plane magneti- 
zation which breaks the in-plane rotational symmetry is 
introduced using the ferromagnetic film, the effect of this 
on the analysis in Sees. MI Al and lHI Bl will depend on the 
angle 9 mentioned above; for instance, a magnetization 
in the y direction will couple to a y and will therefore shift 
the momentum k y for hi and k x for /12. Hence, when ex- 
perimental tests of the various results obtained in those 
two sections are performed, one can probe whether the 
Hamiltonian for the system of interest is actually hi or 
/12 or a linear combination of the two, by varying the 
direction of magnetization of the ferromagnetic film and 
studying the effect that this has on the conductance. For 
example, if the hamiltonian describing the surface elec- 
trons of the topological insulator turns out to be /12, m y 
will have no effect on transport. In general, for any 9, 
there will be specific direction of the in-plane magneti- 
zation m = m y cos 9 + m x sin 9 which will have maximal 
effect on the transport while the component of the mag- 
netization ml — —m x cos 9 + m y sin 9 will not affect the 
transport at all. 
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IV. TRANSPORT IN NMS JUNCTIONS 

We consider a NMS junction on the surface of a topo- 
logical insulator as shown in Fig. [3] As shown there, 
region II, which extends from x — —d to x — 0, has a 
proximate ferromagnetic film leading to an induced mag- 
netization M = My. Region III depicts the supercon- 
ducting region occupying x > d. We assume that super- 
conductivity in this regime is induced via a proximate 
superconducting him with s-wave pairing as shown in 
the figure. The quasiparticles of such a superconductor 
can be described by the following Dirac-Bogoliubov-de 
Gcnncs equation^ 



H — Ep A(r) 
A*(r) Ep — H 



ip = Eip, 



(26) 



where ip = (ipf, ipl, tpV) are the four components for 
the electron and the hole spinors, and the Hamiltonian 
H is given by 



H 



-ihvpa -V + iMJ v B0{x + d)O(-x), (27) 



where 8(x + d) and 9{— x) are Heaviside step functions. 
A(r) = A exp(icj))8(x) is the BCS pair-potential in re- 
gion III. 

Eq. (1) can be solved for the normal, magnetic and 
superconducting regions. In the normal region, the wave 
functions for electron and hole moving in ±x direction 
are given by 

1>f = (l,±e ±lQ ,0,0)exp[z(±fc„x + gy)], (28) 
^ = (0,0,l,Te ±M ')exp[z(±fc^ + gy)], (29) 



sin a = 



Hvpq 



E F V 



sin a = 



twpq 



E F V 



(30) 



where the wave vector k n (k n ) for the electron (hole) wave 
functions are given by 



k n (k n ) — 



Hvf 



(31) 



and a(a ) is the angle of incidence of the electron (hole). 

In region II, the wave functions for an electron and a 
hole moving in the ±x direction are as follows: 



= (l,±e ±i9 ,0,0)eiq>[i{±k b x + qy)], 

= (0,0,l,Te ±i6 ')exp[i(±k' b x + qy)], 
Hv F (q + M) . ,_hv F {q + M) 



sm 



E f 



suit 



(32) 
(33) 
(34) 



where the wave vector k b (k b ) of the electron (hole) wave 
function is given by 



k b (k b ) = 



'e+{-)E F 



-{q + Mf 



(35) 



Here 8(8 ) is the angle of incidence of the electron (hole). 
Note that in principle, we could have applied an addi- 
tional gate voltage Vo in this region as was done in Ref. 
Hil . However, this leads to an expression of the longitu- 
dinal momentum 



k b (k b ) = 




-){E f 



-{q + Mf. (36) 



This shows that in the limit of large Vq, the effect of M 
on k b (k' b ) and hence on G becomes negligible. Therefore 
we restrict ourselves to the Vo = limit. 

In the superconducting region, the BdG quasiparticles 
are mixtures of electron and holes. Hence the wave func- 
tion for BdG quasiparticles moving in ±x directions with 
transverse momenta q and energy e for Ep 3> e, Ao are 
given by 



4? 



sin 7 



( e =Ftf, =Fe ±i(7-/3) ) e -i* Te i(±7- 
x exp[i(±k s x + qy) — kx], 
hv F q 

\E F y ks " 



») 




(37) 
(38) 



and/3 = cos- 1 (e/A )6'(Ao-e)-icosh" 1 (e/Ao)6l(e-Ao) 
where 8 denotes the Heaviside step function. 

Next, we note that for any transmission process to take 
place we need a ,8,8 ,7 < tt/2. This condition gives 
the limits for the range of a. For simplicity we consider 
V = in region II. Then 0' > 6 > a . Using Eqs. (|30]>. 
(f34| and (f38|) . we find that the Andreev process takes 
place for a c \ < a < a C 2, where 

Od = arcsin[(-|e-£; F |)/|e + ^|], (39) 
a c2 = axcsm[(\e- Ep\-M)/\e + Ep\]. (40) 

Note that ot c i 7^ — ck C 2, and this asymmetry is generated 
by the induced magnetization M. 

Following Ref. 1 3, we write wave functions for the nor- 
mal, magnetic and superconducting regions as 



t %j)g, 



(41) 
(42) 
(43) 



where both normal and Andreev reflection are taken into 
account. Here r and ta denote the amplitudes for normal 
and Andreev reflection respectively. These wave func- 
tions must satisfy the following boundary conditions, 



*JvU=-d = *bL=- 



*bL=q = *sU=o- 



(44) 



Solving these boundary conditions, we obtain for 



FIG. 11: Plot of subgap tunneling conductance G(V) as a 
function of the induced magnetization M and applied voltage 
V. See text for details. 



FIG. 12: Plot of zero-bias tunneling conductance G(0)/Go as 
a function of the induced magnetization M and the barrier 
width d. See text for details. 



r,rA,t and t — 



-2ik n d 



N/D, 



(45) 



N = [e ia cos(k b d + 9)+i sm(k b d)} 

-p[cos(k b d -9) + ie ia sm(k b d)\, (46) 
D = [e~ ia cos(k h d + 6)-i sin(fe 6 d)] 

+p[cos{k b d -9)- ie~ ia sin(M)], (47) 
1 



ik n d 



cos(k b d — 



cos(6»)[re-^ + e^] v 
+ie ia sm(k b d)} + re iknd [cos(k b d 
-9)~ie- ta sm(k b d)}), (48) 

t = Tt, (49) 

t'(T+ l)e 2fe " d cos(6>')e-^ 

(50) 



cos(k' b d — 9') — ie ia sin(Avd) 
where the parameters p, T and r\ can be expressed as 



e" 17 - 7? 



e M cos(k b d + 9 ) — i sin(fc 6 rf) 
cos(k' b d — 9') — ie~ la sm(k' b d) 



(51) 
(52) 

(53) 



The tunneling conductance of the NMS junction can 
be expressed in terms of r and ta as 



GjeV) 
G (eV) 



da \l-\r\ 2 ' '-- 12 



cos a \ 

taV cos a. 

cos a / 



(54) 



A plot of the subgap tunneling conductance G/Gq as a 
function of the magnetization M and the applied voltage 
V for a fixed barrier width d is shown in Fig. [TT] We find 
that G(0) decreases monotonically as a function of the 
magnetization for all values of the applied voltage. This 
can be easily attributed to a decrease in the number of 



conduction channels (i.e., number of k y modes with real 
k b ) with increasing M. The behavior of the zero-bias con- 
ductance as a function of the barrier width d and magne- 
tization M is shown in Fig. [T^] We find that the zero-bias 
conductance shows an oscillatory behavior as a function 
of the barrier width d for small M— . With increasing 
M, the position of the conductance maxima shifts which 
demonstrates the tunability of the zero-bias conductance 
with the induced magnetization. This continuous shift 
in position of the zero-bias conductance maxima is to be 
contrasted with the sudden change of its counterpart in 
NMN junctions of topological insulators. 



V. EXPERIMENTS 

Experimental verification of our work would involve 
carrying out the following experiments. For a topologi- 
cal insulator in a crossed electric and magnetic field with 
£ < vfB sin(0), we propose measurement of the energy 
gap of the Landau levels as a function of the electric field 
strength and the tilt angle 9. Such measurements have 
been done in quantum Hall systems using microwave ab- 
sorption techniques 18 . The variation of the excitation en- 
ergy gap between the ground and the first excited states, 
Ai, as shown in Fig. [5] should be observable in similar 
experiments performed with topological insulators. For 
£ > vfB sin(0), we propose measurement of conductance 
G of these films as a function of both the electric field 
£ and the tilt angle 9. We predict that for small £ , the 
tunneling should show a faster suppression with increase 
9 from to ir/2. We note that for a 9 suitably chosen 
between and ir/2, the conductance of these films can 
tuned via an electric field, as demonstrated in Fig. |6l 
leading to realization of electric field controlled switch- 
ing in these materials. For the NMN and NBM junctions, 
which can be prepared by depositing ferromagnetic films 
on the surface of a topological insulator, we propose mea- 
surement of the tunneling conductance G as a function 
of mo. For the geometry shown in the left panel of Fig. 
[21 we predict that depending on the magnetization M, 
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G should demonstrate either a monotonically decreasing 
or an oscillatory behavior as a function of the junction 
width d. Another, probably more experimentally con- 
venient, way to realize this effect would be to measure 
V c of a junction of width d for several values of M and 
confirm that V c varies linearly with M with a slope of 
hvp/(2e), provided /i and d remain fixed. For the geom- 
etry depicted in the right panel of Fig. [21 one would, in 
addition, need to create a barrier by tuning the chemical 
potential of an intermediate thin region of the sample as 
done earlier for graphene^. Here we also propose mea- 
surement of Gi as a function of Vo (or equivalently \) 
for several representative values of mo and a fixed V. 
We predict that the maxima of the tunneling conduc- 
tance would shift from x = n7: to \ = (n + 1/2)tt beyond 
a critical mo for a fixed V, or equivalently, below a crit- 
ical V, for a fixed mo. Finally, for the NMS junction, 
we propose measurement of the tunneling conductance 
G(V) as a function of the magnetization M which should 
demonstrate the decaying behavior shown in Fig.QTJ The 
tunability of the zero-bias conductance maxima, shown 
in Fig. [T2j can also be tested by making junctions with 
different widths. 



In conclusion, we have studied several magnetotrans- 
port properties of Dirac Fermions on the surface of a 
topological insulator, and have shown that they exhibit 
several properties which are distinct both from their 
counterparts in graphene and conventional Schrodinger 
electrons in other 2D systems. These novel features in- 
clude tunability of the orbital and Zeeman effects of an 
applied magnetic field with a crossed in-plane electric 
field, realization of a magnetic switch using a NMN junc- 
tion, and magnetic tunability of transmission resonances 
of Dirac fermions in NBM and NMS junctions. We have 
suggested experiments which can verify our theory. 
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